Universal Quantum Signatures of Chaos in Ballistic Transport 
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The conductance of a ballistic quantum dot (having chaotic classical dynamics and being coupled 
by ballistic point contacts to two electron reservoirs) is computed on the single assumption that 
its scattering matrix is a member of Dyson's circular ensemble. General formulas are obtained 
for the mean and variance of transport properties in the orthogonal (/3 = 1), unitary (f3 = 2), 
and symplectic (/3 = 4) symmetry class. Applications include universal conductance fluctuations, 
weak localization, sub-Poissonian shot noise, and normal-metal-superconductor junctions. The 
complete distribution P(g) of the conductance g is computed for the case that the coupling to the 
reservoirs occurs via two quantum point contacts with a single transmitted channel. The result 
P{g) oc g~ 1+ ^ 2 is qualitatively different in the three symmetry classes. 
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The search for signatures of chaotic behavior in quan- 
tum mechanical systems Q has recently been extended to 
semiconductor nanostructures known as "quantum dots" 
||,||. A quantum dot is essentially a mesoscopic elec- 
tron billiard, consisting of a ballistic cavity connected 
by two small holes to two electron reservoirs. An elec- 
tron which is injected through one of the holes will ei- 
ther return through the same hole, with probability R, 
or be transmitted through the other hole, with probabil- 
ity T. Classically, the uniform (ergodic) exploration of 
the boundaries yields T = R, if the two holes are of the 
same size and sufficiently small that direct transmission 
(without boundary reflections) can be ignored. 

For a closed quantum dot (without holes), it is well 
known that one of the quantum signatures of its clas- 
sically chaotic character consists in the Wigner-Dyson 
distribution of the energy levels [Hj^]. The Wigner- 
Dyson distribution was originally derived by random- 
matrix theory (RMT), and is characterized by a repulsion 
of nearby levels which depends only on the symmetry of 
the Hamiltonian. The quantum analog of the ergodic 
exploration of the dot boundaries by the classical tra- 
jectories consists in the Porter- Thomas distribution of 
the eigenfunctions, as confirmed by numerical studies of 
their amplitude distribution at the boundaries The 
quantum dot with holes is an open, rather than a closed 
system. Just as the Wigner-Dyson distribution describes 
the Hamiltonian H of the closed system, Dyson's circu- 
lar ensemble Q provides the statistical properties of the 
scattering matrix S of the open system. To have spectral 
or scattering properties given by the universal RMT de- 
scription for H or S can be actually regarded as a precise 
definition of the somewhat vague concept of "quantum 
chaos" . To what extent a real ballistic cavity is close to 
this precise universal limit is the subject of the theory of 
quantum billiards |||| . 

Assuming this definition of a quantum chaotic system 



— and this is our only assumption — we will calculate the 
statistics of the transmission and reflection eigenvalues 
of the quantum dot, and hence its transport properties. 
This allows us to determine the universal quantum sig- 
natures of chaos in ballistic transport. Our investigation 
was motivated by a remarkable calculation by P. A. Mello 
of the variance of the conductance in the circular unitary 
ensemble Q . The approach presented below recovers his 
result as a special case, and puts the quantum transport 
theory for a ballistic chaotic billiard on the same footing 
as the established theory for a disordered wire. 

Dyson's circular ensemble characterizes a system where 
all scattering processes are equally probable, subject to 
the constraints of current conservation and time-reversal 
and spin-rotation symmetry. There exist three symmetry 
classes: If a magnetic field B is applied, S is only unitary 
(f3 = 2, unitary ensemble); When B = 0, S is a unitary 
symmetric matrix in the absence of spin-orbit scattering 
(f3 = 1, orthogonal ensemble) or otherwise a unitary self- 
dual quaternion matrix (/3 = 4, symplectic ensemble). 
We assume the usual definition of S, in terms of N x 
N submatrices describing how the incoming modes are 
reflected or transmitted: 

»-(; ?)■ « 

N is the number of transverse modes at the Fermi level in 
each of the two leads connecting the dot to the reservoirs. 
The probability Pp(dS) to find S in a neighborhood dS 
of some given S is 

Pf,{dS) = +r MidS), (2) 

Vp 

where Vp — J /j,p(dS) is the total volume of the S- matrix 
space and fip(dS) is the /3-dependent measure of the 
neighborhood dS of S. In the original work of Dyson , 
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these measures are expressed in eigenvalue-eigenvector 
coordinates. This is a suitable representation to obtain 
the distribution of the scattering phase shifts, but is not 
very convenient for a study of conduction through the 
quantum dot. A transport property A can generally be 
expressed as a linear statistic A = y~] . f(T n ) on the 
transmission eigenvalues T n . The T n 's are not eigen- 
values of S and are not in any simple way related to 
the scattering phase shifts. Instead, T n is an eigenvalue 
of the transmission matrix product tt* . The measures 
fig(dS) have recently been calculated in the transmission- 
eigenvalue representation [0. Since this technical ad- 
vance is at the basis of our analysis, we briefly sketch the 
derivation for the orthogonal ensemble ((3 = 1). 

For P = 1 the scattering matrix is unitary symmetric, 
so that it can be represented in the form S — YY T , 
where Y is unitary. Note that this decomposition is not 
unique. An infinitesimal neighborhood dS of S is given 
by dS — iYdQY T , with dQ a real symmetric matrix. It 
has been shown by Dyson Q , that if the matrix elements 
dQ ij vary through some small intervals of lengths dpij , 
the measure \X\ equals \i\(dS) = Yii<j dfiij, independent 
of Y. We use this freedom to choose Y in the form 



Y = UOT = 







' l + VKU/2 
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where 1Z is an NxN diagonal matrix with elements R n = 
1 - T n . Since idQ = dY T Y* + YUY , and since Y and 
dY can be expressed in terms of the matrices U, O, X 
and their neighborhoods dlA and dO, one can easily get 
dQ in this parametrisation. The result is an expression 
for fii(dS) in terms of the measures n(dU) and fi(dlZ) 
associated with the matrices U and TZ, times a Jacobian: 

Pi(dS) =Y\\Ri-R J \X[{l-R i )- 1 ! 2 Li{dK)ii{dU). (4) 

i<j i 

Integration over the unitary matrix hi gives the reflection 
eigenvalue distribution P(R\, R2, ■ ■ ■ Rn) m the circular 
orthogonal ensemble. The calculations in the unitary fic[ | 
and symplectic pT| ensembles proceed similarly. 

The final result is conveniently written in terms of a 
new set of variables A„ S [0, 00], related to the reflection 
and transmission eigenvalues by R n = A rl /(1 + A„), T n = 
1/(1 + A„). The distribution P(Ai, A 2 , . . . X N ) of the A- 
variables takes the form of a Gibbs distribution, 

P({A n }) = Z- 1 exp[-/3W({A„})], (5a) 

^({A„}) = -£ln |Ai W*)' ( 5b ) 
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l<3 

2/3 



ln(l + A), 



(5c) 



where Z is a normalization constant. The symmetry pa- 
rameter (3 £ {1,2,4} plays the role of an inverse tem- 
perature. The fictitious "Hamiltonian" H consists of a 



logarithmic pairwise interaction plus a one-body poten- 
tial V/3(A). This potential is symmetry-independent to 
order N, while the term of order A^ depends on (3. 

Remarkably, the distribution (|J) is identical to the 
global maximum-entropy ansatz for the transfer matrix 
of a diffusive conductor jL2j — except for the one-body 
potential, which is different: The potential Vd(A) for a 
disordered wire of length L and mean free path I is given 

by pl 



V d (A) = (Nl/L)ln 2 (V\ + VTTX) +O(N ). 



(6) 



The potential (J5c[), in contrast, contains no microscopic 
parameters and increases more slowly with A. In the case 
of a disordered wire, it is known |H| that the logarith- 
mic repulsion — In | A; — Xj \ is only rigorously valid for the 
weakly reflected scattering channels (Aj, A^ <C 1). In the 
ballistic chaotic dot, the logarithmic repulsion which we 
have found is a direct consequence of the basic assump- 
tion that the scattering matrix belongs to the circular 
ensemble. It remains to be seen whether this assump- 
tion breaks down in some range of A. It seems that 
this problem is unrelated to the breakdown of the log- 
arithmic repulsion in the statistics of the energy levels of 
a chaotic system JIJJl^], which occurs on energy scales 
greater than the Thouless energy E c . This energy scale 
has no obvious counterpart for the transmission eigenval- 
ues. 

We consider transport properties of the form A = 
Sn=i a (^«)- To calculate the expectation value (A) = 
L a(X)p(X)d\ for a ballistic chaotic system, we need the 
density p(X) of the A's in the circular ensemble. For this 
purpose, we use Dyson's large- TV expansion JTJ] 



A - A' 



dX' 



where -f denotes the principal value. We decompose 
P = Pn + Sp into a contribution pat of order N (giving the 
"Boltzmann conductance" ) and a symmetry-dependent 
correction 5p of order N° (responsible for the "weak- 
localization effect"). From Eqs. (j^) and ([?]) we find, order 
by order, 1 



dX = r => Pn(X) 



N 



A- A' 



1 + A 



7r(l + A)VA' 



(8) 



1 The result (g) for the O(N ) correction Sp holds only 
for A <C N 2 / 3 , because for larger A's the O(N) contribu- 
tion pM ~ A r A _3/ ' 2 no longer dominates the density and 
the large-A^ expansion fails. The large-A tail ensures that 
J °° Sp(\) dX = 0, but is irrelevant for the conductance of the 
quantum dot. [The range A > N 2 ^ 3 3> 1 gives a contribution 
to g of order N~ 2 ^ 3 , which can be neglected relative to the 
contribution of order 1 which is retained.! 
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Sp(X) 



0-2 
4/3 



MA), (9) 



where the one-sided delta-function satisfies 
/ °°* + (A)dA = l. 

It is interesting to compare this eigenvalue density for a 
ballistic chaotic cavity to the density p^ for a disordered 
wire. To order N and for L ^> I one has [fLTfl 
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2L ^A(l + A) 



, for A < A c ~ ±exp(2£/Z). 



(10) 



The density goes to zero abruptly near a cutoff A c S> 1, 
in such a way that pff(X) dX — N. The transmis- 
sion eigenvalue density pn(T) — pN(X)\dX/dT\ (with 
T = (1 + A) -1 ) is bimodal, with a peak near unit and 
near zero transmission. This is a familiar result for a dif- 
fusive conductor, but was not previously established for 
a chaotic dot. 

The term of order A^ , which yields the weak- 
localization (anti-localization) corrections for a disor- 
dered wire when /3 ^ 2, is given by [|l8fl: 



V d) (A) = 



ff-2 
2/3 



<5+(A) + (A + A 2 )- 1 / 2 



x (41n 2 [VA + VTTA]+tt 2 ) '1 , (11) 



which is not as strongly peaked near A = as the delta- 
function result ([)]) for a chaotic dot. 

We now use Eqs. (@) and (0) to calculate the expec- 
tation value (T) of the total transmission probability 
T = Trttt = £) n (l + A„) _1 . According to the Landauer 
formula, T equals the conductance g (measured in units 
of 2e 2 /h where the factor 2 comes from spin or Kramers 
degeneracy). The result is 



(T) = ±N + 6T, ST=^ 



(12) 



For (3 = 2, one finds (T) = \N = (R) (where (R) = 
N — (T) is the total reflection probability). This is the 
quantum analog of what we expect from the "ergodic" 
exploration of the dot boundaries by the classical trajec- 
tories. Quantum interference then breaks the equality 
(T) = (R) by an amount ST, due to weak localization 
(/3 = 1) or anti-localization (/3 = 4). The value ST = — j 
for /3 = 1 is in agreement with the results of Ref. (l9) and 
demonstrates the point raised in Ref. (2(J : Weak localiza- 
tion is not only given by coherent backscattering, but has 
an off-diagonal (in mode index and classical trajectory la- 
bels) component. In the same way one can compute the 
average of any other linear statistic. We give two exam- 
ples to illustrate the generality of our approach. 

The first example is the shot-noise power P which 
is given by P = P Trift(l - with P Q = 
(2e 2 /h)2ell (U is the applied voltage). In this case 



a(A) = P A(1 + A)" 2 . Since Sp(X)a(X) = for any /3, 
there is no weak-localization correction for the shot noise 
of a chaotic dot — in contrast to a diffusive conductor, 
where a weak- localization effect does exist [E2| . Integra- 
tion of pn(X)o,(X) gives the average shot-noise power 



(P) = lNP = iPpoisson, 
o 4 



(13) 



which is four times smaller than the Poisson noise 
-Ppoisson = gPo — 2el associated with a current / of un- 
corrected electrons. The \ reduction in a chaotic dot is 
to be compared with the ^ reduction of shot noise in a 
diffusive conductor |23f| . 

The second example is the conductance Gns of the dot 
if one of the two attached reservoirs is a superconductor. 
This case corresponds to @ a(A) = (4e 2 //i)(l + 2A)~ 2 
if /3 = 1,4 (Gns is not a linear statistic for (3 = 2). 
Again, we find a noticeable difference between the dif- 
fusive disordered wire and the ballistic chaotic dot. In 
the disordered case, the conductance Gn in the normal 
state is unchanged if one of the reservoirs becomes su- 
perconducting ((Gn) = (Gns)) u P to a weak-localization 
correction of order N°). In the ballistic chaotic case, Eq. 
(|) yields (to order N): 



(G NS ) = ^-(2-V2)ZV, 
h 



(14) 



which differs to order A^ from the result (Gn) = 
(2e 2 /h)(N/2) in the normal state. Another example, 
which we do not discuss here but which can be treated 
in just the same way ]2^| , is the Josephson effect if both 
reservoirs are superconducting. 

So far we have focused on the expectation values in 
the circular ensemble. Fluctuations around the average 
in this ensemble can be computed using the general for- 
mulas of Ref. pi| , which hold for any ensemble with a 
logarithmic interaction (regardless of the form of the one- 
body potential). The variance in the large- A^ limit is 
given by 



VarA = -— — 



1 1 

pV 2 



dX 



da(A') 



x In 



VA - y/X' 



(15) 



Eq. (15) can be used to compute the analogue of the 
"Universal Conductance Fluctuations" (UCF) in a bal- 
listic chaotic cavity. These fluctuations can be induced 
by a change of the Fermi energy, of the applied magnetic 
field, or of the boundary spin-orbit scattering, as well as 
by a slight deformation of the shape of the dot. One 
obtains, for example, Varg = 1/8/3, Var P/Po — 1/64/3 
for the fluctuations in the conductance and shot noise, 
respectively. The l//3-dependence is the same as for a 
disordered wire, but the numerical coefficients are some- 
what different due to the difference in interaction poten- 
tial ||. 
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The results for mean and variance given above require 
iV > 1. The opposite regime N — 1 is also of in- 
terest. This would apply to a semiconductor quantum 
dot which is coupled to the reservoirs by two quantum 
point contacts with a quantized conductance of 2e 2 /h. 
The probability distribution (||) reduces for N = 1 to 
P(A) = \(3{l + A)~ 1_/3 / 2 . This implies for the (dimen- 
sionless) conductance g = (1 + A) -1 the distribution 



P(g) = \Pg- 1+0/ \ o <<?<!. 



(16) 



This is a remarkable result: In the presence of magnetic 
field (/3 = 2), any value of the conductance between 
and 2e 2 /h is equally probable. In non-zero held it is 
more probable to find a small than a large conductance, 
provided that the boundary scattering preserves spin- 
rotation symmetry {(3 — 1). In the presence of spin-orbit 
scattering at the boundary (/? = 4), however, a large con- 
ductance is more probable than a small one. To observe 
this qualitatively different behavior presents a challenge 
for experimentalists. 

In summary, we have calculated the distribution of 
the transmission and reflection eigenvalues characteriz- 
ing Dyson's circular ensemble. Relying on a definition of 
"quantum chaos" based on the applicability of this en- 
semble to describe scattering in a ballistic chaotic cavity, 
we have extracted from the joint probability distribu- 
tion (||) the expectation values and the fluctuations of 
arbitrary linear statistics related to quantum transport. 
We have mainly stressed the differences between ballistic 
chaotic dots and disordered wires. Another important 
point, which we underline in conclusion, consists in the 
qualitative differences existing between a fully chaotic dot 
and a dot where the classical dynamics is integrable and 
for which larger quantum fluctuations are generally ex- 
pected [g§. 
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Note added: Upon completion of this manuscript we 
received a preprint by Baranger and Mello, in which some 
of our results are obtained by a different method. 
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